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Abstract 

We show the existence of an infinite family of finite-time singularities 
in isotropically expanding universes which obey the weak, strong, and 
dominant energy conditions. We show what new type of energy condition 
is needed to exclude them ab initio. We also determine the conditions 
under which finite-time future singularities can arise in a wide class of 
anisotropic cosmological models. New types of finite-time singularity are 
possible which are characterised by divergences in the time-rate of change 
of the anisotropic-pressure tensor. We investigate the conditions for the 
formation of finite-time singularities in a Bianchi type VIIo universe with 
anisotropic pressures and construct specific examples of anisotropic sud- 
den singularities in these universes. 

PACS numbers: 98.80.Cq, 98.80.Bp, 98.80.Jk 



1 Introduction 

It has recently been shown ^ |2] that the Einstein equations permit the de- 
velopment of finite-time singularities during the evolution of an expanding uni- 
verse that obeys the strong energy condition. In isotropic and homogeneous 
cosmologies this singularity is characterised by a divergence in the pressure, 
p, and the acceleration of the expansion scale factor a as t — > t<; < cxd. The 
density p, expansion scale factor a, and expansion rate a/ a all remain finite 
at ts- In order to exclude this type of singularity it is sufficient to introduce 
a pressure-boundedness condition, for example that p < Cp for some finite 
positive constant C, or to require that dp/dp be continuous ^ Without 
such a condition linking the pressure to the density, a pressure singularity can 
occur independently of the density and leads to a violation of the dominant 
energy condition (expressed by \p\ < p) as t ^ tg ^E]. Note that the diver- 
gence of p with finite p leads to an unbounded quantity with dimensions of a 
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velocity. Similar forms of singularity can be obtained in Friedmann universes 
arising as cosmological solutions of higher-order gravity theories |2| E| and other 
forms of finite-time singularity have been proposed in the presence of cosmo- 
logical accelerations of fluid flow relative to the hypersurface-orthogonal fluid 
congruence in |^. Sudden singularities are sp curvature singularities that are 
neither strong-curvature nor crushing [2] and their modest effects on geodesies 
have been studied in ff^. These mathematical features of simple cosmological 
models impinge upon investigations of the late-time behaviour of the so called 
'phantom matter' models. The latter generally have unconventional behaviour 
leading to future 'big rip' singularities at flnite time more severe than those dis- 
cussed here jS]. It should be noted, however, that not all phantom cosmologies 
end in a big rip singularity 0. 

We recall that, in the case of isotropic universes, sudden singular behaviour 
results from a scale factor of the form ^ 

a(t)=(l)\«.-l) + l-(l-l)", (1) 

with Us = a(ts). Hence, as t ^ ts from below, we have 

a^,i,-l)Bt^-^-^j^0^j^^-oo, (2) 

whenever 1 < n < 2 and < g < 1. This solution of the Friedmann equations 
exists in the interval < t < ts- More generally, a finite-time singularity will 
arise in a Friedmann universe when [2] 

a{t)^(^.Ly ^as-l) + lHt,-tri^pJ^^^ , (3) 

where ajk are constants, Nj < j are positive integers, Q £ Q"*" and the quantity 
in {..} brackets is a convergent double psi series. The latter tends to zero as 

t^ts. 

In this paper we extend previous work on the subject by studying a new 
family of sudden singularities in isotropic Friedmann universes and the forma- 
tion of sudden singularities in the presence of anisotropic expansion. In ^U] an 
anisotropic and inhomogeneous cosmology of Stephani type was found to pos- 
sess similar finite-time singularities, as well as spatial analogues that preserve 
the dominant energy condition within finite regions of space. We shall investi- 
gate whether new forms of anisotropic finite-time singularity can arise, possibly 
under weaker energy conditions than in the isotropic case. 

We begin our discussion with a covariant formulation of the sudden singu- 
larity problem in section 2. This allows us to identify the various possibilities 
that may occur. In section 3 we describe a new infinite family of sudden sin- 
gularities in isotropic universes which obey all the standard energy conditions. 
In section 4 we give an explicit example of some of the possibilities uncovered 
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in section 2 by constructing anisotropic sudden singularities in a Bianchi type 
VIIq universe with anisotropic pressures. We summarise with our conclusions 
and a discussion of finite-time singularities in section 5. 



2 Covariant formulation 

We begin by investigating under what conditions new types of finite-time singu- 
larity can arise in an anisotropic, irrotational, spatially homogeneous, spatially- 
flat spacetimes, containing matter with density p, isotropic pressure p and a 
trace- free anisotropic pressure tensor iTab- The evolution of this model is de- 
scribed by the following set of covariant equations (e.g. see [lip 

p = -e(p + p)-a„b^''^ (4) 

e -ie^ - i(p + 3p) -2cr2, (5) 

CTah = — fQcTafc - <J c{a(^'^ b) ~ Eab + ^T^ab , (6) 
Eab = -QEab — ^{P + P)crab - ^T^ab - ■^'dT^ab + Sac/^aE'^b) 

-^CTcia-^'b) ■ (7) 

Here, is the volume expansion scalar, aab is the shear tensor, describing kine- 
matical anisotropics, and Eab is the electric part of the Weyl tensor, which is 
associated with tidal forces. We proceed by assuming that p is independent of 
p and that iTab is independent of both p and p. We will also leave the time evo- 
lution of ITab unspecified by any constitutive relationships. These propagation 
equations must be supplemented by the constraint equations 

92 = 3{p + <J^), (8) 

O'ab = — ©CTafc + TTab, (9) 

both imposed by the spatial flatness of the model. Note that the first of these is 
the generalised Friedmann equation and the latter is the residual Gauss-Codacci 
formula, which here provides an alternative description of the shear evolution. 
Also, expressions © and © combine to give 

TTab = ^QcTab - 2ac{aCr'^b) " ^Eab , (10) 

which directly relates all the sources of anisotropy. 

Introducing non-zero spatial curvature leaves Eqs. Q)-® the same, although 
generally it adds a cmlHab term to the right-hand side of (UJ (where Hab is the 
magnetic Weyl tensor). In the presence of curvature Eqs. lO and © also change 
into 

7^ = 2(p-ie2 + o•2) , (11) 

and 

T^{ab) = -<^ab ~ (^(^ab + T^ab , (12) 
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respectively. Here, TZ(ab) is the symmetric and trace-free component of TZab, the 
spatial Ricci tensor, TZ = TZ°'a is the associated Ricci scalar and cr^ = craf,(T°^/2 
is the magnitude of the shear tensor. Also, Eq. (|12|l combines with Q to give 

T^{ab) = -jQcTab + <Jc{a'^'^b) + ^ab + ^T^ab ■ (13) 

According to the system (B|)- (|1U|I . we can have an anisotropic- pressure sin- 
gularity where the quantities a, O, p, p and Uab are all finite but -Kab diverges 
as t tg < oo. Here the singular behaviour is manifested as a finite time 
singularity in (Tab, the first time derivative of the shear tensor (see Eqs. ® 
and (O).^ Let us now look at this possibility from a different perspective. 

Consider the scale factor a, defined by means of the isotropic volume ex- 
pansion, according to a /a = 0/3. Then, Raychaudhuri's formula (see Eq. jSJ) 
transforms into 

''^ = -l{p + 'ip)-la\ (14) 
Applied to a Bianchi / spacetime, the time derivative of the above gives 

= -iep - ip + iQcjabcr"' - laabT^"" , (15) 

given the spatial flatness of the model. In an anisotropic cosmology with non- 
zero spatial curvature, however, this generalises to 

^ = -l^P -\P+ fe^ab^"' " hab^""' + |^a67^<'^^> , (16) 

where TZt^ab) satisfies expression Following 1)14 11 - 1)16(1 . we see that a diverges 
when TTab is not bounded. At the same time the lower-order derivatives of the 
scale factor remain finite. It worth noticing that, unless appropriate regularity 
conditions are imposed on the matter, a diverging p will also lead to a singularity 
in a (see Eq. (fTTHl V 

When dealing with anisotropic expansion it helps to consider the scale factor 
defined by means of the generalised Hubble law. In covariant terms the latter 
reads J^l 

^ = ^ + a„fce"e^ (17) 

where Co are unitary, linearly independent, spacelike vectors (i.e. Cae" — 1 and 
Caw" = 0). The shear term in the right-hand side of the above ensures that, in 
contrast with a, the new scale factor also contains information on the anisotropy 
of the expansion. The full role of the anisotropy is revealed by taking the time 
derivative of (|17|l . On using (0 the latter leads to 

^ = ^ + (cTabe-e") ' - + E^b - i^afc) e'^e^ , (18) 



^One may also consider the vacuum case by setting p = = p = iTab- Then, the reduced 
system of Eqs. IH- llijl appears to allow for a finite time singularity in a^f, if 'R-i^ab) or ^ab 
diverge. However, when is bounded within the time interval in question, a series of known 
existence theorems due to Wald and Rendall prevents this from happening )12| . 
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with aja satisfying expression ifT^ . In addition, one can use expression lfT7|) to 
define the individual scale factors along the three linearly-independent directions 
determined by (with i = 1, 2, 3). For example, contracting Oah twice along ei 
in ll7|l we obtain 

which defines the scale factor in the ei-direction. Similarly, Eq. H18|l leads to 



- = - + cr 

ii a 
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- En + i^ii , (20) 



with analogous expressions for £2 and £3. Results (|18|) and H20() show that, when 
""afaj Eab or TZ(^ab) diverges, £ and £i will generally diverge as well, which is in 
distinct contrast with the behaviour of a. Recall that an unbounded anisotropic 
pressure leads to the singular behaviour of a' instead of a (see Eqs. (|14ll . (|15|l ). 
Note also that, unlike a diverging nab, a diverging p manifests itself as a singu- 
larity in the third time derivative of £ and not in the second. 

So far, our inspection of the self-consistent divergences in the sets (jH)- ()10|) . 
(O-ldni and |(ni)-(j2nil has provided necessary rather than sufficient conditions 
for these new types of finite-time singularities to occur. Clearly, their detailed 
form must also comply with all of the Einstein field equations. In what follows 
we will construct explicit examples to show how some of the aforementioned 
situations may arise. 



3 New isotropic examples 

It is clear from Eq. (|15|) that in the isotropic case (aab — 0) there exists an 
infinite family of Fricdmann sudden singularities with the scale factor a(t) chosen 
to have the form or its generalisation Q, which satisfy the weak energy 
condition (p > and p + p > 0), the strong energy condition {p + 3p > and 
p + p > 0), and the dominant energy condition (p > and p ± p > 0). If we 
choose n in Eq. ^ so that 

N < n < N + 1 , 

where iV is a positive integer, then there will be a finite-time singularity in the 
{N + ly* time derivative of the scale factor. In other words, 

fitN+i ^ °° as t^ts, 
with all the lower order derivatives remaining finite, namely 

-^j^ L < 00 as t ^ts, 
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for all r = 1, 2, . . . N. In this limit the (A^ + 1)** time derivative of the scale 
factor produces an infinity in the {N — 1)** time derivative of the pressure at 
ts, but the density p and all lower derivatives oi p remain finite as i — > is since 

1 d^+ig _i d^-y 

Specifically, for a scale factor of the form we see that 



V ^-N^^ 14-U MW+l-ii ' ^ 



If we choose iV = 1, we create the sudden singularities given in |2] and 
also by Eq. (|2Jl above. Because they arise in d and -p they lead to a violation 
of the dominant energy condition despite satisfying the weak and strong energy 
conditions. However, the infinite family of sudden singularities characterised 
by the choice N >2 satisfy the weak, strong and dominant energy conditions 
and produce no divergent spacelike energy fluxes of the type discussed in 
Unbounded behaviour occurs in the time derivatives of the pressure and might 
suggest the need to define a family of new energy conditions if this behaviour is 
to be excluded by fiat. For example, we could exclude finite-time singularities of 
this type if we introduced a generalised matter regularity condition that required 

d^ ^ ^ , . 

- — < Cs- — , for some s < r , 

dr At^ - 

where Cs is a positive constant. 



4 Anisotropic examples 

Let us look in more detail at some of the possibilities revealed in the last section. 
Consider the spatially homogeneous, anisotropically expanding, Bianchi VIIq 
universe with a line element given by .14. J. 5., , 16, 

ds^ = ~At^ + g„^dx"dx^ , (21) 

where t is the comoving proper time and 

(A^ [cosh /x-h sinh /xcos(fc2;)] sinh/xsin(fcz) \ 

sinh/xsin(fcz) [cosh /x — sinh /icos(fcz)] (22) 

j 

with a, /3 = 1, 2, 3 and z = x^; also, A = A{t), B — B{t), fi — ii{t) and fc is a 
constant parameter. The above metric can be interpreted as the superposition 
of an axisymmetric Bianchi I spacetime and a circularly polarised gravitational 
wave propagating along the axis of symmetry (see j^El and also QB]). The 
wavenumber and amplitude of this wave is given by k and /z respectively. Note 
that the above metric has anisotropic spatial curvature as well as anisotropic 
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expansion and reduces to the axisymmetric Bianchi / spacetime in the long- 
wavelength limit (i.e. as /c — > 0). It is the most general spatially homogeneous 
generalisation (or perturbation) of the zero-curvature Friedmann universe and 
has played an important role in discussions of the isotropy of the universe |17| . 

When matter is in the form of a fluid with anisotropic pressures, the Ein- 
stein equations associated with the metric (|21f) . (|22|l lead to the evolution for- 
mulae [13 CS]: 



A A A B\ ^. 



(23) 



B 



AB 



o+2To=^(P + P3-pi-P2) + ihj sinhV 



and 




's) sinh(2^) = _pi - P2 ■ 



These are supplemented by the following constraint: 




— ] sinh^ II 
B ' ^ 



= P- 



(24) 
(25) 

(26) 



which is the generalised Friedmann equation. The above reduces to the standard 
Friedmann equation for the flat FRW metric, where A — B and /i = 0. In 
Eq. (|26l) A = A{t) and B = B{t) represent the two scale factors that characterise 
the anisotropic expansion. Also recall that ^ = ^(t) and k are respectively 
the amplitude and the wavenumber of the superimposed gravitational wave of 
wavelength 2TTB/k. 

Consider now the following power-law evolution equations for the two indi- 
vidual scale factors 



A{t) = As-1 
B{t) = Bs-l 



- A{ts - tr , 



+ Bits-ty., 



(27) 
(28) 



where ts corresponds to a time in the late evolution of the model, Ag = A{ts), 
Bs = B{ts) and A, B are constants. Also, g, m G M and, for reasons that will 
become clear next, we will assume that 2 < n and that 1 < r < 2. Then, taking 
the first and second time derivatives of A we find 



A = 



ts \t 



9-1 



+ nA{ts - ty 



A = ^(^("IV '-n(n-l)^(t,-t)"-^ 



(29) 
(30) 
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while by successively differentiating B we obtain 

/ , \ m— 1 

B = ^^fc^(i)"" + r(.-W.-0^-. (32) 

Given that 2 < n, Eqs. H29|) and (|3U|I guarantee that all of A, A, and A remain 
finite as t approaches ts- On the other hand, since 1 < r < 2, expressions (I31f) . 
H32|l imply that B diverges as t ts, while B and B remain finite. In addition, 
since £ = {AB'^Y^^ is the average scale factor of the anisotropic expansion, both 
£ and £ are finite at the t ^ tg limit but £ ^ oo. All these properties ensure 
that at t ~ tg the model described by Eqs. lf77|) and experiences a sudden 
singularity at tg of a sort discussed in section 2. Assuming the metric ()21f) . 
(|22|) . this sudden singularity is triggered by a divergent anisotropic pressure at 
tg. For example, when one or both of the principal pressures pi and p2 diverge 
as t ^ tg, expression (|24|1 guarantees that B also diverges at tg. Note that the 
singular behaviour of 13 is in agreement with the power-law evolution of H28|l. At 
the same time, Eq. (|23|) guarantees that A remains finite, which is compatible 
with definition H27|) . Also, according to (|25|l . the divergence of either pi or p2 
immediately implies a singular behaviour for jl, unless pi — P2- Clearly, the rest 
of the variables, namely p, P3, fi and /i, do not show any singular behaviour and 
all the terms in Eq. 1)26(1 are always finite at i > 0. Since p remains finite at 
tg, while at least one of the principal pressures diverges there, leading also to a 
divergence of p, we conclude that in our example the dominant energy condition 
is violated at the t = tg limit ^ 0] . 

Alternatively, we may assume that both pi and p2 are finite, and allow for 
Ps to diverge as t tg. Then, Eqs. H23|l and 1)24(1 ensure that both A and B 
diverge at tg, while 1(25(1 shows that fl remains bounded at that point. This type 
of singular behaviour is also described by the scale factors ((27(1 and 1(28(1 . with 
the extra proviso that now n must also lie within the open interval (1,2). 

It is straightforward to generalise these examples in the same way that we 
did with isotropic singularities in the last section. By choosing n and r to lie 
in appropriate open intervals bounded by successive positive integers we can 
produce mutatis mutandis an infinite family of sudden singularities that occur 
in arbitrarily high derivatives of one or both of the anisotropic scale factors and 
in the associated time derivatives of nab- 

Finally, we note that the anisotropic scale factors given in Eqs. ((27(1 - 1(28(1 
can be generalised in the same way that the particular isotropic solution (^ was 
generalised using the logarithmic series to expression 10 in ^ . More specifically. 
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we have 



^^^^ [i=Ofc=0 

B(t) - 1 + + B{t, -tylY^Yl ^^^(^^ - ty^'i^'^s'its ~t])\ , 



j=0 k=0 



where m,n,q, and r are as before and Ujk, bjk are constants, Nj,Mj < j are 
positive integers and Q,S & Q+. 



5 Conclusions 

We have shown how a new infinite family of finite-time singularities can arise 
in isotropically expanding universes which obey the strong, weak and dominant 
energy conditions. This is possible because they produce infinities in the time 
derivatives of the pressure of arbitrarily high order. In particular, singularities in 
d^+^a/dt^+^ arise from singularities in d^~-^p/dt^"-^ for all positive integers, 
N. These can only be excluded by introducing new regularity conditions on the 
derivatives of the pressure. 

We have also shown that finite-time singularities can arise in new ways when 
the expansion of the universe is anisotropic. The simple forms of sudden singular 
behaviour found in isotropic universes are not rendered unstable by expansion 
anisotropics. Moreover, we have found that new types of finite-time singularity 
can arise from the anisotropic part (nab) of the fluid pressure. We investigated 
the role of anisotropic pressure in detail. Beginning with a general discussion 
of the issue of anisotropic finite-time singularities, we identified a number of 
possible singular cases. We constructed a specific anisotropic model of Bianchi 
type VIIq to show how an unbounded anisotropic pressure tensor can lead to 
finite-time singularities. The latter can occur whilst the matter density and 
isotropic pressure remain finite and can be avoided only by introducing appro- 
priate restrictions on nab- For example, one might require that TTafoTr"'' < Cp^ 
for some positive constant C. In summary, our investigations have shown that 
finite-time singularities are a feature of a wide class of cosmological solutions to 
the Einstein equations and are not simply a consequence of special symmetries. 
In order to exclude them from studies of the general behaviour of the Einstein 
equations at late times we need to introduce new restrictions on the allowed 
form of the energy-momentum tensor. 
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